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We examine phase separation in aqueous mixtures at fixed amounts of hydrophilic monovalent 
ions. When water is the minority component, preferential solvation can stabilize water domains 
enriched with ions. This ion-induced precipitation occurs in wide ranges of the temperature and 
the average composition where the solvent would be in one-phase states without ions. The volume 
fraction of such water domains is decreased to zero as the interaction parameter \ (dependent on 
the temperature) is decreased toward a critical value for each average composition. 
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In fluid mixtures containing water, phase separation 
behavior can be drastically changed by a small amount 
of a salt More strikingly, many groups have ob- 

served long-lived heterogeneities (sometimes extending 
over a few micrometers) in one-phase states and a 
third phase visible as a thin plate at a liquid-liquid inter- 
face in two-phase states Very recently, mesophases 
with lamellar or onion structures have been found for an 
antagonistic salt composed of hydrophilic and hydropho- 
bic ions [3| • Dramatic ion effects are ubiquitous in various 
soft matters. For example, in polyelectrolytes, the phase 
behavior sensitively depend on the degree of ionization 
and the composition for mixture solvents [5[ and large- 
scale heterogeneities have been observed [6]. In these 
phenomena, the solvation interaction among ions and 
polar molecules should play a major role together with 
the Coulomb interaction among charges Q- Recently 
some theoretical efforts have been made to elucidate the 
solvation effects in phase transitions in electrolytes and 
polyelectrolytes |8l-ll0|. In this Letter, we consider hy- 
drophilic monovalent ion pairs such as Na + and Cl _ in a 
binary mixture of water and a less polar component (oil 
or alcohol) and examine ion-induced precipitation. 

Neglecting the electrostatic interaction but accounting 
for the solvation interaction, we first consider a binary 
mixture in a cell with a fixed volume V. Here ions consti- 
tute the third component with density n(r). The volume 
fractions of water, oil, and ions are written as (p(r) and 
(p'(r), and Ujn(r), respectively, where uj is the ionic vol- 
ume. If the two solvent species have the same molecular 
volume i>o, their densities are (p/vo and (p' /vq. The space- 
filling condition reads (p + <f>' + Vin = 1. The ion volume 
fraction is assumed to be small or vju <C 1, which is eas- 
ily satisfied for small ions with vi <C vq. In this Letter 
we thus set <f>' — 1 — <f>, which simplifies the calculations. 

When the ions have a strong preference of water over 
oil, we set up the free energy as 



f = / <lr 



T 



n ln(m>o) — 9 n( P 



(1) 



where C is a positive constant and g is a positive param- 
eter representing the ion preference of water. The space 
integral is within the cell. The Boltzmann constant will 



be set equal to unity and then the temperature T repre- 
sents the thermal energy of a particle. The free energy 
density f((p) is taken to be the simple form, 



v f/T = c(>\n<j>+(l 



iln(l 



+ X#1 ~4>), (2) 



where x is the interaction parameter dependent on T and 
its mean- field critical value is 2 without ions. We fix the 
total particle numbers of the three components as 



drn/V, <f> = / drcj)/V, 



(3) 



In equilibrium the homogeneity of the ion chemical po- 
tential 8F/8n yields 

n = ne 90 /< e90 ). ( 4 ) 

where (e 9 ^) = J dre 9 ^/V is the space average of e 9 ^. 
Substitution of Eq.(4) into F in Eq.(l) gives 



F 
T 



dr 



T 



Vn\n[nv Q /{e a % (5) 



In equilibrium the chemical potential difference h - 
5F / 54> for the composition is also homogeneous. Here, 



h = f'((j>) - TCV 2 4> - Tgn 



(6) 



where /' = df/d(f>. Around a planar interface varying 
along the z axis, we obtain TC{d<f)/ dz) 2 = 2H(cp), where 

H(<t>) = f{(j>) - f(<t>*) - T(n - n a ) - h{<P - a ). (7) 

We suppose coexistence of a water-rich phase a and an 
oil-rich phase j3 with A<f) = <fi a — 4>p > 0. The composi- 
tions and the ion densities in the two phases are written 
as <p a , (pp, n a , and rip, respectively. The volume fraction 
of the phase a is denoted by 7 Q . Since that of the phase 
/3 is jp — 1 — 7 a , we have 



7o 



(pp)/A(p = {n- rip)/ An, 



(8) 



where An = n a — rip. From Eq.(4) the ratio of the bulk 
ion densities is written as n a /np — exp(gAcp) 3> 1 from 



2 



Eq.(4), where we assume gA(j> 3> 1. Neglecting the sur- 
face free energy, we express F as 



F 
V 



7a/c* + 7£//3-Tnlog[( 7a e^ + 7pe a ^)/hv ], (9) 



where f a = f(<fi a ) and fp — f{<j)p). We minimize this F 
with respect to j a , <p a , and <pp at fixed to obtain 

h = fa ~ T 9 n a = fp- Tgnp, (10) 
f a -fp- TAn = hA^, (11) 

where f a — f(<fr a ) and fa — f{(f>p). These equations 
also follow from Eqs.(6) and (7). 
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FIG. 1: Compositions <f) a and (f)p (left) and y a , ry a n a /n, and 
n a /n (right) vs x, where n — 6 x lO" 4 ^ 1 and g = 11. For 
4> = 0.35 (top), 0c continuously changes for \ > Xp — 0.687. 
For </> = 0.46 (bottom), 4> a jumps at \ — 2, where \p — 1-834. 
Shown also are the coexistence curve without ions and the 
spinodal curve with ions. The latter follows from a shift of 
the spinodal curve without ions by g 2 von/2 to the left [f|. 



Our main finding is that there appears a precipitation 
branch of <f> = <j) a = 1 in the region \ < 2 for j > 1. In 
Eq.(10) we set f a = v^ 1 T[- log(l - <j> a ) - x] to obtain 

l-0a = A Q exp[- 9 G(^)], (12) 

where Aq, = exp(x + vq/UT). Here we introduce 

Gf» - -«o[/(0) + (1 - </>)f'(<f>)]/T 

= -lo g( /»-x(l-^) 2 , (13) 

where the second line follows from Eq.(2). By setting 
f a = in Eq.(ll) we also find v n a — G(<pp). Out- 
side the spinodal curve without ions, we notice dG/d(j) = 
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FIG. 2: x = Xp(0) for v n = 4x 10~ 4 , 2x 10" 4 , and 5x 10~ 5 . 
Precipitation occurs for cj> < cj> p (arrows). Shown also are the 
coexistence and spinodal curves without ions. 



-v (l - (j})d 2 f /d<j) 2 < 0, which leads to G(0) > from 
G(l) = 0. Thus the right hand side of Eq.(12) is small 
or <j) a = 1 for large g. The equation for <pp is given by 

Wo n Q = G(<f>p) = ^ o nA0/[0 - ^ + e'^Acf)] (14) 

where we may set A(j) = 1 — (f>p. As j a — ((j) — (f>p)/A(f> 
decreases to zero, (f>p and n a increase up to <f> and 
vone 9 ^ 1- ^, respectively. Remakably, 4> a and <pp depend 
on <f>. From the second line of Eq.(13), the precipitation 
branch exists only for \ > Xp (</>)■ The lower bound is 



Xp 



log(l 



i;one 9(1 -^]/(l - 4>) 2 



(15) 



where 7 Q ->■ as \ -> Xp- 

To easily understand the mathematics, let us focus on 
the case "f a <C 1, where f a = and fp = /(0) — — 
4>)"f a in Eq.(9). Treating j a as an order parameter, we 
write the deviation AF = F(j a ) — F(Q) of the free energy 
from the one-phase value F(0) as 



AF/TV 

-i. 



Ai7 Q -nlog(l + Bu a ), 



(16) 



where A 1 = VqG{4>) + ng(l - (j>) and B x = e 9 ^ - 1 S 
e 9A ^. For to = hBi/Ai > 1, Ai 71 has a negative minimum 
given by — TVAi(w\ogw — w + l)/B\ < attained at 



la 



(w-l)/B 1 ^v n/G($)-e- 9 ^\ (17) 



which is consistent with Eq.(14) in the limit <pp — > cj). 
The condition w > 1 is equivalent to x > Xp f° r 5^1- 
Figure 1 gives the phase diagrams in the <f>-x plane 



with n — 6 x 10 



Sr 1 



and g = 11. In the first case of 



4> = 0.35, 4> a changes continuously and is minimum at 
X = 2.05, where Xp — 0.687 and the maximum of n a is 

X P - 



o.asiw^ 1 at x 



In the second case of 



0.46, 



where Xp = 1-834 and n a = 0.0842WQ - at x = Xp ; we 
find that <f> a changes discontinuously along a hysteresis 
loop in the range 1.953 < x < 1.967. In equilibrium, F 
is minimized and the resultant discontinuous transition 



3 



is at x = 1-957. In Fig. 2, we display curves of x — Xp(0) 
for three values of n with g = 10. Each curve assumes a 
minimum at small <fi far away from the coexistence curve 
without ions. With increasing 0, the precipitation branch 
shrinks and disappears as <j> — > <f> p , where <p v is a critical 
composition less than 0.5. As functions of <f> and Xj we 
show 7 a in Fig. 3 in the continuous case of g = 10 and <p a 
in Fig. 4 in the discontinuous case of g = 11. 

Without the electrostatic interaction, the surface ten- 
sion of our system is expressed as the integral a — 
J dzCTCdtj) ■/ 'dz) 2 around an interface varying along the 
z axis @- Use of Eq.(7) gives 

a = (2CT) 1 / 2 d^H^) 1 / 2 , (18) 

J 9 ft 

where C is assumed to be a constant. In Fig. 5, we display 
the function [2iJ(0)w o /r] 1/2 for g = 11. Here H(<f>) k 
/"O^OO - </>/3)72 as -> 0^ with /" = d 2 f/d<j> 2 . Thus 
we obtain cr - TC(A0) 2 /2£, where £ = (f'/CT)- 1 / 2 is 
the correlation length at = 0^ = 0. 




FIG. 3: Volume fraction of the water-rich phase 7 a as a func- 
tion of x and 0. for n = 2 x lO -4 ^ 1 and g — 10. In this case 
7 a decreases continuously to zero as \ Xp- 




FIG. 4: (Color on line) Composition of the water-rich phase 
<p a for vofi = 6 x 10 -4 and g = 11. In this case a discontinuous 
transition occurs for x — 1-95 and > 0.395 (in red). One- 
phase states are realized for X < Xp or f° r 4> > 0p — 0.473 
and x < 2 (in white). For x > 2, Q is little affected by ions. 
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FIG. 5:_ [2H((f>)v /T} 1/2 vs 0, vanishing at = <f> a and tj> p , 
where — 0.46, n = 6 x lO -4 ^ 1 , and g — 11. 



Including the electrostatic interaction, we next con- 
sider aqueous mixtures containing a hydrophilic mono- 
valent salt. The cation and anion densities are written 
as ni and ri2, whose total amounts are fixed as 



J drni = J drn-2 = Vn/2. 



(19) 



The electric potential $ satisfies the Poisson equation 
V • eV$ = — Aive(ni — 712), where the dielectric constant 
e(0) can depend on 0. The free energy F reads Q 



F 
T 



dr 



f{4>) , C 



|V0| : 



e|V$| 



T 2 8ttT 



(20) 



The ion chemical potentials due to solvation, written as 
/i* ol (</>), strongly depend on 0. The interaction terms 



(oc gi) in F follow for the linear forms /i* ol (</>) = [i^ — Tg^ 
(where the first terms are irrelevant constants) . This lin- 
ear dependence is adopted to gain the physical conse- 
quences in the simplest manner. For each ion species i, 
the solvation-chemical-potential difference between the 
two phases is given by A/i^ = T<j , jA0, which is the 
Gibbs transfer free energy in electrochemistry In 
aqueous solutions, gi ^> 1 for hydrophilic small ions, 
while gi < for hydrophobic ions [1, [ll[ . We minimize 
F with respect to n$ under Eq.(19) to obtain 



11, 



n 



°exp[^Te$/T], 



(21) 



,0 



where — is for i = 1, + is for i = 2, and n\ 
n/[2 J drexp(gicf) =p e$/T)]. The composition profile is 
determined by the homogeneity of SF/S<j>. Here we ne- 
glect the image interaction, whose role is reduced com- 
pared to that of the solvation interaction for not small 
ion densities [1, [l2| . 

For gi = gi — g, we have n\ = = n/2 and <f> = 0, 
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FIG. 6: Left: Numerical droplet radius d/a vs x (bold line) 
together with the theoretical curve without the surface free 
energy (dotted line). Right: Normalized potential e<E>(r)/T 
and water volume fraction <j>(r) (top), and normalized ion 
densities voni(r) and vori2(r) (bottom), where \ — 1-7. Here 
gi = 15, gi — 7, <j> = 0.35, and won = 6 x 1CF 4 . 



so F in Eq.(20) reduces to F in Eq.(l). For gi ^ g 2 . 
an electric double layer appears at the interface with a 
potential difference A$ = T(gi — .Q2)A0/2e across it, but 
the bulk phase relations (10) and (11) still hold with 

<?=(<?! + <? 2 )/2. (22) 

In Fig. 6, we give numerical results for g\ — 15, gi = 7, 
<f> — 0.35, and vqu = 6 x 10~ 4 . The space unit is a — v^ 3 . 
We suppose a spherical water-rich droplet with radius d 
placed at the center of a spherical cell with radius R — 
600a. Then 7 Q = (d/i?) 3 . The dielectric constant is of 
the form e = 40(1 + 4>). We also set aC = x and e 2 /T — 
120a. In the left, the droplet disappears at d = 26.0a, 
where w in Eq.(17) is 1.08. This critical radius follows 
if the bulk free energy AF in Eq.(9) is equated with the 
minums of the surface free energy 4irad 2 . In the right, we 



set x = 1.7 to display the profiles of ^(r), </>(r), ni(r), and 
n2{r). We obtain (f> a = 0.993 and n a = 0.352v ~ 1 within 
the droplet and 4>p = 0.349 and rip — 2.55 x lO -4 ^ 1 
outside it. In Fig. 6, the potential $(r) relaxes with the 
Debye length k^ 1 — 11.8a and is well fitted to the one- 
dimensional solution of the nonlinear Poisson-Boltzmann 
equation (dotted line) with A$ = 3T/e Q. Here n a and 
np are the bulk values of n = n\ + ri2- 

In future we should explain the experimental findings 
of large-scale heterogeneities Q • We note that one-phase 
states are metastable outside the spinodal curve with ions 
(X < 1/20(1 — (f>) —g 2 vofl/2) in Fig. 1. Thus precipitation 
from a one-phase state should be triggered by some im- 
purities and/or hydrophilic walls. We also note that the 
wetting transition of aqueous mixtures is much influenced 
by the ion-induced precipitation mechanism. 

Experiments are informative, where the temperature, 
the water volume fraction, and the salt amount are var- 
ied. We mention an experiment by Leunissen et al. 
fl3T |. where micron-sized water droplets containing ions 
formed a crystal in an oil with low dielectric constant 
(e ii = 4 — 10) without a surfactant. Graaf et al. 
ascribed its origin to the screened Coulomb interaction 
among droplets. We also propose experiments of the 
salting-out effect of polyelectrolytes in water-alcohol Q, 
where the degree of ionization much increases with accu- 
mulation of water around the polymers [lOj . 
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